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Abstract 

The approximate stress-energy tensor of the conformally invariant massless spin- 1/2 field in the 
Hartle-Hawking state in the Schwarzschild spacetime is constructed. It is shown that by solving the 
conservation equation in conformal space and utilizing the regularity conditions in a physical metric 
one obtains the stress-energy tensor that is in a good agreement with the numerical calculations. 
The back reaction of the quantized field upon the spacetime metric is briefly discussed. 
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Recently, the renormalized stress-energy tensor of the quantized conformally invariant 
massless spin- 1/2 field in the Schwarzschild and Reissner-Nordstrom spacetimes has been 
evaluated numerically . This important and long awaited result completes our knowl- 
edge of the behaviour of three basic quantum test fields in the Hart le- Hawking state. How- 
ever, obtained as a by- product of the calculations its analytical approximation is, contrary 
to the spin and spin 1 cases, very poor near the event horizon (although it is accurate at 
large distances). 

It has been shown 0, 0] that the stress-energy tensor of the massless fields in the 
Schwarzschild geometry decomposes naturally as 

{T£) rm = (T^) analytic + A£, (1) 

where (T^) anal y tlc is the analytical approximation of the stress-energy tensor and A^ is a 
traceless and covariantly conserved tensor that must be calculated numerically. For the 
conformally invariant massless scalar field the analytical part coincides with the Page ap- 
proximation (and the approximations constructed in (], 0, [J), whereas for the vector 
fields such approximation has been constructed in Refs. |^, |8( . It has been shown that for 
both scalar and vector fields the approximation is reasonable. (Analytical approximation of 
tensor Af, of the scalar field has been constructed in kjl). 

, (T^) analytic for the spin- 1/2 field gives 
a wrong sign of the energy density at the event horizon, invalidating thus any prospect ap- 
plications. On the other hand, however, the expectation value of the stress-energy tensor of 
the conformally invariant massless fields in the Hartle- Hawking state in the Schwarzschild 
geometry is known to possess some general features [10]. The asymptotic behavior of tan- 
gential and radial components of (TJf) ren , the regularity conditions on the event horizon and 
the trace anomaly are quite restrictive, and allow construction of a class of approximate 
tensors. Further, a piece of the numerical data, such as the exact value of one of the com- 
ponents of the stress-energy tensor, say (Tg) ren , on the event horizon may be used in the 
final determination of the model. It should be noted that calculations of the horizon value 
of the stress- energy tensor could be regarded as a relatively simple task since it is sufficient 
to retain only the two smallest frequencies in the mode sums 

an 

In this sense this 

information could be regarded as independent of the numerical calculations of the stress- 
energy tensor for r > 2M. 
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The idea of reconstructing (T^) ren from the knowledge of its asymptotic behaviour is not 

I 

new and bel ong s to Christensen and Fulling [10]. It has been subsequently elaborated in 



Refs. 



In this note we shall show that, contrary to the widespread opinion, it is 



possible to construct the stress-energy tensor of massless spin 1/2 field in the Schwarzschild 
rich satisfactorily approximates the 'exact' {Tjj) ren . The method is similar to 



spacetime, w' 

that of Ref 16| with the one reservation: here we consider the four component spinors 
rather than the two component ones. Moreover, we shall explicitly demonstrate that the 
Frolov-Zel'nikov method |?| also yields reasonable results and subsequently compare the 
stress-energy tensors constructed within the frameworks of both methods. 

Let us start by counting available information. First, it is known that the stress-energy 
tensor is covariantly conserved and its trace is given by a general formula 

(T p ) ren = a (H + OR) + /3Q + 7 QR, (2) 

where 

H = R p , vpT R iiVpT — c lR ia/ R! lv + g-R 2 ; (3) 

g = R^ pr R^ pT - AR^nr + r 2 (4) 

and the numerical coefficients for the spin-1/2 field are a — 18 A, (3 — —11 A, 7 = 0, and 
A = 2 7 457r 2 . Further, we observe {T p ) ren is regular on the past and future event horizon 
and approaches at large distances the flat spacetime radiation stress-energy tensor. Finally, 
we expect that the curvature effects enter the stress-energy tensor as ~ x 3 (x = 2M/r). 
This requirement is usually motivated by the observation that far from the event horizon 
the stress-energy tensor should consist of the red- shifted thermal bath part supplemented 
by the quantum corrections. This assumption is sometimes referred to as a weak thermal 
bath hypothesis [17| . 

The idea is to construct T{? in the optical metric and subsequently to transform it back to 
the physical spacetime with the aid of a transformation that relates the stress- energy tensor 
in conformally related geometries 0, 0, ^| . Now, the independent informations listed above 
suggest that the tangential component of the stress-energy tensor could be approximated as 



^l + f>a^ 



f e = o T 1 + , (5) 



with iV = 6 and T" 1 = 90vr 2 (8M) 4 . 

Returning to the Schwarzschild geometry and making use of the regularity conditions in 
the physical metric, after some algebra, one has 



9 2 



1 + 2x + 3x 2 + —x 3 + -(230 + Ua,)x A + —x 5 + —x 6 
7 14 v ' 7 7 



(6) 



Remaining components of the approximate stress-energy tensor can be easily obtained from 
© and will not be displayed here. Similar results for the two component spinors have been 
constructed in Ref. [lS] to which the reader is referred for the technical details. It should be 
noted that the logarithmic term x 2 (ai + 202) Inx, which appears as the result of integration 
of the conservation equation survives even if the regularity conditions are satisfied. Only 
after accepting the weak thermal bath hypothesis the coefficients a\ and a 2 could be equated 
to zero |8(. The parameter 04 can be determined form the equation: 

r;(2M) = (t;(2M)>, (7) 

which, by the spherical symmetry, equality of T\ and T r r at the event horizon and the trace 
anomaly yields in fact only one condition. 

Although the horizon values of the components of the stress-energy tensor in the Hartle- 
Hawking state have never been cited explicitly in literature, and Carlson et al. present 
their results only graphically, we have sufficient informations to construct the model. The 
horizon value of the trace anomaly in the case at hand is 7/76807r 2 M 4 , whereas the approx- 
imate values of T[f is jjj 

T e J2M) = T£(2M) w 1 -r . (8) 

eK ' ^ ' 90tt 2 (8M) 4 V ; 

Making use of Eqs. (JHJ) and Eq. (JJJ) with the right hand side given by (JSJ) one finally obtains 
T\ = -Ir^3 + 6x + 9x 2 + 12x 3 -^ + ^^-69x 6 ^ , (9) 

7 „ n 2 52s 3 139a; 4 18s 5 15x 6 \ 
K = -T[l + 2x + ?>x 2 - — + — — + (10) 

. 7 m / o 68x 3 89x 4 102x 5 87x 6 \ , , 

T e e = -Til + 2x + 3x 2 + — _ + + (11) 



and 



The run of the components of are displayed in Figs. 1-3. The new 'radial' coordinate £ 
is defined as £ = (r — r + )/M, where r + denotes location of the event horizon. 



It should be emphasized that the case of the spin-1/2 fields is special. Indeed, for the 
electromagnetic field the T" coincides with the analytic part of (T") |8j, whereas the analo- 
gous result constructed for the conformally invariant scalar fields substantially improves the 
Page approximation. 

A different method of calculating the stress-energy tensor in static spacetimes has been 
proposed by Frolov and Zel'nikov in Ref. Q|. It has been shown that it is possible to 
construct a family of expressions describing the approximate stress-energy tensor T^ Z ' V 
solely form the curvature, the Killing vector and their covariant derivatives up to some 
given order. By construction the Frolov- Zel'nikov tensors have a correct trace and a proper 
behaviour under the scale transformations. Upon imposing appropriate regularity conditions 
at the event horizon of the Schwarzschild black hole one obtains a one parameter family of 
the approximate stress- energy tensor. As expected, the resulting expressions are simple 
polynomials in x. If for the conformally invariant scalar fields the free parameter is set to 
zero the result coincides with the Page approximation. It should be noted, however, that 
there are no a priori reasons to accept such a choice. Indeed, for a vector field such a 
choice leads to evidently wrong results. On the other hand the free parameter may be easily 
adjusted employing Eq. (J2J). 

It could be easily shown that 

T (FZ)v = T u + ^ (12) 

where is the approximate stress-energy tensor (jHHUJ) and the conserved and traceless 
tensor D u is given by 

t _ 7 /783X* 174x 5 g7x 6x 



7 /232x 3 667 x 4 174a; 5 261 x 6 . 

D r = -T\ 1 1 (14) 

2 1 35 35 35 35 ' 



and 



» 7^ /116a; 3 145a; 4 174a; 5 87 x 6 , 

Dl = — T + + . 15 

2 1 35 7 35 7 1 



The components of D^ are displayed in Figs. 1-3 (right panel). 

Now we can compare the approximate tensors (piTTj) to (|12H15|) and to the numerical 
results presented in Ref. Q]. By construction the tensors are exact at infinity and very close to 
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the exact value at the event horizon. Inspection of similar figures presented in Ref. |2j shows 
that our approximation (fillip is in good agreement with the exact numerical calculations. 
The Frolov-Zel'nikov approximation is slightly worse but still reasonable. Unfortunately, 
we are unable to provide detailed comparison between the exact numerical results and the 
approximations as the former were presented only graphically. 

The run of and T^ F constructed for the massless spinor field qualitatively resembles 
behaviour of the approximate stress-energy tensor of the electromagnetic field 4J. Indeed, 
for the electromagnetic field the Frolov-Zel'nikov approximation is also less accurate than 
the analytic part of (T^). 




FIG. 1: This graph shows the radial dependence of the rescaled component T\ [A = 90(8M) 4 7r 2 ] of 
the approximate stress-energy tensor of the massless spin 1/2 field in Schwarzschild spactime and 
the rescaled component D\ (right panel). 

Now we shall employ the constructed approximation in the backreaction calculations. It 
should be noted however, that in order to investigate the evolution of the system consisting 
of a black hole and the quantized field it is necessary to have at one's disposal the detailed 
knowledge of the functional dependence of the renormalized stress-energy tensor on a wide 
class of metrics. So long the effects of the quantum gravity could be safely ignored this is 
achieved solving the semi-classical equations: 



(16) 



The first term in the right hand side of the above equation comes from the classical source 
whereas the second one is due to the contribution of the quantized fields. Unfortunately, 




FIG. 2: This graph shows the radial dependence of the rescaled component T£ [A = 90(8M) 4 7r 2 ] of 
the approximate stress-energy tensor of the massless spin 1/2 field in Schwarzschild spactime and 
the rescaled component D r r (right panel). 




1 2 3 4 5 



FIG. 3: This graph shows the radial dependence of the rescaled component Tjj [A = 90(8M) 4 7r 2 ] of 
the approximate stress-energy tensor of the massless spin 1/2 field in Schwarzschild spactime and 
the rescaled component D@ (right panel). 

all our current understanding of the stress-energy tensor of the quantized massless fields is 
limited to the static black holes described by the electric charge and the mass. Our results, 



however, may be of use in the linearized back reaction calculations [19|, |20j: the curved 
spacetime leads to the stress-energy tensor approximated by (JSHHJ), which, in turn, modifies 
the background geometry. 

Since the stress-energy tensor is asymptotically constant it is necessary to put the system 



7 



in the cavity or spherical box of definite radius, r c . It should be noted, that presence of 
the boundary certainly modifies the stress-energy tensor. Therefore the radius r c should 
be chosen to guarantee applicability of the perturbative approach on the one hand and to 
minimize the error caused by ignoring expected boundary effects on the other. 

The geometry of the quantum corrected spherically-symmetric static black hole is gener- 
ally described by the line element 

ds 2 = _ e 2^r) f ^ dt 2 + J_ dr 2 + r 2 dn 2^ 

where f(r) = 1 — 2m(r)/r and dQ 2 is the metric on a unit sphere. We shall assume that the 
functions m(r) and ijj(r) can be expanded as 

m(r) = M + eMi(r) + 0(e 2 ) (18) 

and 

V>(r) = eVi(r) + 0(e 2 ). (19) 

To keep control of the order of terms we have introduced the dimensionless parameter e, 
which will be set to 1 at the final stage of calculations. Similarly, (T"[g]) in the right hand 
side of Eq.(Jinj) should be substituted by e{T"[g]). 

The solutions of the linearized semi-classical Einstein field equations with a source term 
given by the stress-energy tensor (p iTTj) reduce to two elementary quadratures: 

Mi(r) = -4tt J r 2 T l t dr + C x (20) 

and 

Mr) =4,/ + C, (21) 

where the integration constants C\ and C2 are to be determined form the boundary condi- 
tions. Observe that it is possible to determine the function ipi(r) without prior knowledge 
of the horizon value of the stress-energy tensor. It is simply because the difference T r r — T/ 
does not depend on the coefficient a^. 

Our preferred choice of the boundary condition for Eq. ()20|) is simply Mi(r + ) = 0, which 
requires knowledge of the exact location of the event horizon, r + . From ()18|) and (|20|) one has 
r + = 2M, and hence M is to be interpreted as the horizon defined mass. For the function 
■01 (r) we shall adopt the natural condition gtt(r c )g rr (r c ) = —1. 
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Now the equations (|2Uj) and (}2~T]) can be easily integrated to yield 



K 

M 1 (r) = — 



7 



7 



21 



13 



31 



M V6x 3 
161 , 



+ 



6 



-x 



2x 2 ' 2x 
14 lnx 



33 H x x 



(22) 



and 



K 

M 2 



1(1 

6 V x 2 



+ 7 



x 



25 



-— [x 2 - x 2 c ) - 13 ( x 3 - xfj - 14 In 



1 

X c 
X 
X (■ 



50 



(x — x r 



(23) 



where K = 38407T and x c = 2M/r c . It is believed that the line element ()17|) with the 
metric potentials given by Eqs. (jl8pi9j) with (|2^jl and (|2~3j) respectively better approximates 
the physical reality then the original (Schwarzschild) one. Having the quantum corrected 
geometry of the black hole one may study its properties such as Hawking temperature, trace 
anomaly and its influence on the motion of test particles. Since the calculations of these 
effects are elementary we shall not dwell on them here. We only remark that in view of the 
results of Ref. Q] presented method is the only one, which is able to provide simple and 
reasonable approximations to the exact stress-energy tensor. This model may be thought of 
as a minimal one, i. e. requiring only one numerical information: the horizon value of the 
stress-energy tensor, which, as has been pointed out, is to certain extend independent of the 
numerical calculations for r > 2M. Of course, more complicated models which approximate 
the exact {Tp better could be easily devised at the expense of the additional numerical 
informations. 

Similar calculations with different asymptotics may be used in construction of the stress- 
energy tensor of the spin-1/2 field in the Unruh state. We expect that the approximation will 
be reasonable and the results presented in Refs. suggest that it is a safe anticipation. 
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